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' Abstract. Our numerical simulations of the Edgeworth-Kuiper belt ob- 

jects gravitationally scattered by the four giant planets accounting for 
mean motion resonances reveal numerous resonant gaps in the distribu- 



q . tion of the scattered population. 



\ Available numerical simulations indicate that gravitational scattering of the 

EKBOs by the four giant planets might explain the transport of comets from 
the trans-Neptunian region all the way inward, down to Jupiter, as well as the 
origin of the so called scattered disk (for a review, see Malhotra et al. 1999). Our 
recent simulations exploring both the spatial and phase structure of cometary 
populations beyond Jupiter (Ozernoy, Gorkavyi, & Taidakova 2000 = OGT) 
deal with 36 test bodies, initial positions of which were taken close to the known 
EKBOs so that their orbits intersected the Neptune's orbit. Positions of every of 
the scattered test bodies were traced for ~ 0.5 Gyrs and recorded each revolution 
of Neptune about the Sun, e.g. 165 yrs, totalling ~ 2.8 • 10 7 positions. This 
delineates in great detail the spatial as well as phase structure of the cometary 
population beyond Jupiter. 

Our simulations reveal that (i) each giant planet dynamically controls a 
cometary population called the "cometary belt" (a well-known 'cometary fam- 
ily', such as the Jupiter family comets, is just a visible part of the cometary belt); 
and (ii) comets avoid resonant orbits so that the belts contain numerous gaps in 
the (a, e)- and (a, i)-spaces similar to the Kirkwood gaps in the main asteroidal 
belt. We have found that a numerous resonant gap structure is mantained in 
the scattered cometary population far beyond the Neptunian belt. 

Fig. 1 shows how the distribution of the scattered bodies in semimajor axis, 
a, evolves on a time scale of 0.5 Gyrs. Every bin of Aa = 0.1 AU was visited by 
each of the test bodies, on average, ~ 500 times, which confirms a strong level of 
chaos in the scattered population. Our following findings are worth mentioning: 
(i) Numerous gaps at the resonances with Neptune are well pronounced; 
(ii) A robust resonant structure of gaps is rapidly established; (iii) Unlike to the 
distribution of the scattered objects in heliocentric distance, which has a maxi- 
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Figure 1. Evolution of the distribution of simulated scattered 
comets in semimajor axis: a for comets with 25 < a < 50 AU and 
b for comets with 50 < a < 100 AU. Number of comets is given be- 
tween a and a + Aa, where Aa = 0.1 AU. Arrows indicate various 
resonances with Neptune. Ten distributions, from below to the top, 
are shown at 1, 2, 4, 8, 16, 32, 64, 128, 256, and 500 Myrs, respectively. 



mum near the Neptunian orbit (Levison & Duncan 1997, OGT), the distribution 
in semimajor axis appears rather uniform; (iv) There are appreciable accumu- 
lations of the scattered comets near the resonances well outside the zone of a 
strong gravitational scattering (e.g., near the 2:1 resonance). Those accumula- 
tions (called 'diffusive' ones in OGT) are produced by slowing down of diffusion 
as the comet's distance of pericenter turns out to be outside the planetary orbit. 
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